Unconventional isotope effects in the high-temperature cuprate superconductors 
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We review various isotope effects in the high-T c cuprate superconductors to assess the role of the 
electron-phonon interaction in the basic physics of these materials. Of particular interest are the 
unconventional isotope effects on the supercarrier mass, on the charge-stripe formation temperature, 
on the pseudogap formation temperature, on the EPR linewidth, on the spin-glass freezing temper- 
ature, and on the antiferromagnetic ordering temperature. The observed unconventional isotope 
effects strongly suggest that lattice vibrations play an important role in the microscopic pairing 
mechanism of high-temperature superconductivity. 
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I. INTRODUCTION 

Developing a correct microscopic theory for high-T c su- 
perconductivity is one of the most challenging problems 
in condensed matter physics. More than ten years af- 
ter the discovery of the high-T c cuprate superconductors 
by Bednorz and Muller [jj], there have been no micro- 
scopic theories that can describe the physics of high-T c 
superconductors completely and unambiguously. Due to 
the high T c values and the earlier observation of a small 
oxygen-isotope effect in a 90 K cuprate superconduc- 
tor YBa2Cu3C>7_y many theorists believe that the 
electron-phonon interaction cannot be the origin of high- 
Tc superconductivity. Most physicists have thus turned 
their minds towards an alternative pairing interaction of 
purely electronic origin (e.g., see Ref. §-§). As a matter 
of fact, the idea that the highest T c is only 30 K within 
the conventional phonon-mediated mechanism is not jus- 
tified. In principle, T c can increase monotonically with 
both the phonon frequency and the electron-phonon cou- 
pling constant within the phonon-mediated Eliashberg 
theory ||. The limitation of the T c within the phonon 
mechanism is only imposed by a possible structural insta- 
bility in the case of too strong an electron-phonon inter- 
action. However, there is no universally accepted, simple, 
and quantitative stability criterion ||. 

It is well known that the observation of a gap in the 
electronic excitation spectrum || and the discovery of 
an isotope effect [tTojPl in conventional superconductors 
provided important and crucial clues to the understand- 
ing of the microscopic mechanism of superconductivity. 
In particular, the effect of changing isotope mass on the 
superconducting transition temperature T c implies that 
superconductivity is not of purely electronic origin, but 
that lattice vibrations (phonons) play an important role 
in the microscopic mechanism for this phenomenon. 

The first evidence for an isotope effect was reported 
in 1950 by Maxwell JTy] and independently by Reynolds 
et al. pT| . They found that the critical temperature T c 
of mercury is an inverse function of the isotope mass. 
In the same year Frohlich |12] pointed out that the same 



electron-lattice interaction which describes the scattering 
of conduction electrons by lattice vibrations gives rise to 
an indirect interaction between electrons. He proposed 
that this indirect interaction is responsible for supercon- 
ductivity. Frohlich's theory got strong support from the 
observed isotope effect, and played a decisive role in es- 
tablishing a correct mechanism. In 1956, Cooper fTs] ] 
demonstrated that electrons with an attractive interac- 
tion form bound pairs (so called Cooper pairs) which lead 
to superconductivity. However, the existence of electron 
pairs does not necessarily imply a phonon mediated pair- 
ing. Indeed, Bose condensation as considered in 1955 by 
Schafroth [Q, is also a possible mechanism for super- 
conductivity, but the model was not able to explain the 
isotope effect. Finally, in 1957, Bardeen, Cooper and 
Schrieffer developed the BCS theory which can ex- 
plain most physical properties observed in conventional 
superconductors. 

Remarkably, the BCS theory can well explain the iso- 
tope effect. The T c within the theory is given by 



fcflT c = 1.137kJDexp 



N(0)V 



(1) 



where wjj is the Debye frequency, which is proportional 
to M^ 1 / 2 . The electron-phonon coupling N(0)V is the 
product of an electron-phonon interaction strength V and 
the electronic density of states at the Fermi surface N(0), 
both of which are independent of the ion mass M in the 
harmonic approximation. Eq. 1 implies an isotope-mass 
dependence of T c , with an isotope-effect exponent a = 
—dh\T c /d\nM = 1/2. This is in excellent agreement 
with the reported isotope effects in the non-transition 
metal superconductors (e.g., Hg, Sn and Pb). In fact, the 
isotope effect was the first justification of the proposed 
electron-phonon coupling mechanism. 

The conventional phonon-mediated superconducting 
theory is based on the Migdal adiabatic approximation 
in which the phonon-induced electron self-energy is given 
correctly to the order of (m^/M) 1 / 2 ~ 10~ 2 , where mb 
is the bare mass of an electron. Within this approxima- 
tion, the density of states at the Fermi level N(0), the 
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electron-phonon coupling constant A ep , and the effective 
mass of the supercarriers are all independent of the ion 
mass M. However, if the interactions between electrons 
and nuclear ions are strong enough for electrons to form 
polarons (quasiparticles dressed by lattice distortions), 
their effective mass m* will depend on M. This is be- 
cause the polaron mass m* — rrif, exp(A/ ui) Jl6| , where 
A is a constant, and w is a characteristic optical phonon 
frequency which depends on the masses of ions. Hence, 
there is a large isotope effect on the carrier mass in po- 
laronic metals, in contrast to the zero isotope effect in 
ordinary metals. The total exponent of the isotope effect 
on m* is defined as j3 = J2 -d\mn* / dh\Mi (Mi is the 
mass of the ith atom in a unit cell). From this definition 
and the expression for the polaron mass m* mentioned 
above, one readily finds 



ln(m* /rrib). 



(2) 



The above equation implies that there should be a large 
negative isotope effect on m* in polaronic metals. 
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Therefore, if the electron-phonon interaction is strong 
enough to form polarons and/or bipolarons, one will ex- 
pect a substantial isotope effect on the effective mass of 
carriers. In this article, we will present various isotope 
effects in cuprates including the unconventional isotope 
effects on the supercarrier mass, on the charge-stripe for- 



mation temperature, on the pseudogap formation tem- 
perature, on the EPR linewidth, on the spin-glass freez- 
ing temperature, and on the antiferromagnetic order- 
ing temperature. These unconventional isotope effects 
clearly demonstrate that phonons are relevant to the ba- 
sic physics of cuprates and may be important for the 
occurrence of high-temperature superconductivity. 



II. ISOTOPE EFFECT ON THE 
SUPERCONDUCTING TRANSITION 
TEMPERATURE 

Studies of isotope shifts of T c have been carried out 
in almost all known cuprates. A comprehensive review 
was given by Franck |l7| ] . The role of the anharmonicity 
of the apical oxygen on the isotope effect was discussed 
in detail by Miiller |H| . Most of the studies reported so 
far concern the oxygen-isotope shift (OIS) of T c by re- 
placing 16 O with 18 O, partly because the experimental 
procedures are simple and reliable. Now it is generally 
accepted that optimally-doped cuprates exhibit a small 
and positive oxygen- isotope exponent ao- It was also 
found that ao for optimally-doped materials decreases 
with increasing T c [|Tg|] . The small OIS observed in the 
optimally doped cuprates suggest that phonons might 
not be important in bringing about high temperature su- 
perconductivity. 
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FIG. 2. The oxygen 
exponents as a function of T c for Yi_ :E Pr a ;Ba2Cu307-i / 
Yi_ a; Pr 2 :Ba2Cu408. The data were taken from li"" 

However, the doping dependence of the OIS has 
been extensively studied in different cuprate systems 
]l7| , p0| -|2l| . For a particular family of doped cuprates the 
OIS increases with decreasing T c , and can be even larger 
than the BCS value. In Fig. |l|a, we plot the doping depen- 
dence of T c for the 16 and 18 samples of the single-layer 
cuprate I^-zSr^CuO^ It is clear that the T c 's of 18 
samples are always lower than those of the 16 samples. 
The doping dependence of the isotope exponent ao is 
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shown in Fig. |T]b. The magnitude of ao increases with a 
decrease of doping and becomes very large (> 0.5) in the 
deeply underdoped regime. The large ao value observed 
near x — 0.125 might be related to the structural insta- 
bility [p8[ . The results suggest that the phonon modes 
related to the oxygen vibrations are strongly coupled to 
conduction electrons. 
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FIG. 3. The site-selective oxygen-isotope effect in 
YBa 2 Cu 3 06.96. Here < 18 O p 18 O ac > means all the oxygen 
sites are exchanged by ls O, and < 18 O p 16 O ac > means that 
only the in-plane oxygen sites are replaced by ls O. It is ev- 
ident that the planar oxygen mainly (> 80%) contributes to 
the total oxygen-isotope shift in the optimally doped cupratc. 
After lES. 



In addition to the large oxygen-isotope shifts observed 
in underdoped compounds, there are also large copper 
isotope shifts observed in underdoped La2_ x Sr a .Cu04 
Hi, oxygen-depleted YBa 2 Cu 3 7 _ y |||, Pr-substituted 
Y 1 _ x Pr x Ba 2 Cu 3 7 _j / |§ and Y 1 _,Pr :c Ba 2 Cu 4 8 p0|, 
as well as in YBa2Cu40s 



31 1 . In Fig. 0, we plot 



the oxygen and copper isotope exponents as a function 
of T c in Yi_ x Pr x Ba2Cu307_ y and Yi_ 2; Pr x Ba2Cu408. 
As T c or doping decreases, both ao and acu increase 
monotonously. Interestingly, ac u is about 3/4 of ao in 
the deeply underdoped region, while ac u is even larger 
than ao near the optimal doping. This suggests that 
the Cu-dominated phonon modes (most of them are low- 
energy modes) are involved in the superconducting pair- 
ing. 
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FIG. 4. The site-selective oxygen-isotope effect in an opti- 
mally-doped YBa 2 Cu 3 6 .94. After Q. 

The large copper-isotope shift also implies that the 
phonons in the CuC>2 planes are relevant to supercon- 
ductivity. By analogy, one should also expect that the 
planar oxygen vibrations make more important contri- 
butions to the pairing than the apical and/or chain oxy- 
gen vibrations. We can distinguish the contributions 
of the different oxygen sites to the total OIS by site- 
selective oxygen-isotope experiments. A partially site- 
selective oxygen-isotope exchange was attempted by Car- 
dona et al. J|^] and by Ham et at. |33|. The experimen- 
tal difficulty in doing a complete site-selective oxygen- 
isotope exchange was overcome by Nickel et al. f54j. In 
their experiment the authors replaced the 18 O by 16 O 
in the chain and apical oxygen sites of a fully 18 O ex- 
changed YBa 2 Cu3C>7 sample, while keeping the 18 un- 
exchanged in the plane sites. The site-selectivity was 
confirmed by Raman spectroscopy Q . A small negative 
OIS [ao = —0.010(4)] associated with the planar oxygen 
was found. In contrast, Zech et al. fl35| showed a small 
positive OIS [ao = 0.018(4)] related to the planar oxy- 
gen ions. The above discrepancy may be due to a broad 
superconducting transition in the samples of Nickel et al. 
p4j , which makes it difficult to define T c reliably. 

Although the OIS is small in the optimally-doped sam- 
ples, the experiment done by Zech et al. |35]| provided 
evidence that the planar oxygen ions mainly (80%) con- 
tribute to the total OIS (see Fig. ||). Nevertheless it is 
not obvious that the same conclusion should apply to 
the underdoped samples where the total OIS is large. 
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Zhao et al. |36|] thus carried out the site-selective oxygen- 
isotope experiments in the underdoped and optimally 
doped samples of Yi_ K Pra;Ba2Cu307- a . In Fig. f|, we 
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phonon coupling? It is very unlikely that the isotope 
effect is due to a difference in the carrier densities of two 
isotope samples. This is because the T c 's of the 18 sam- 
ples are always lower than the 16 O samples by more than 
1 K in the underdoped region, independent of whether 
dT c /dx is positive, negative, or zero (see Fig. |l]a). 
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FIG. 6. The total (T) isotope shifts as a function of T c for 
Yi_ a! Pr :E Ba2Cu307_ y together with the isotope shifts from 
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chain (AC) oxygen ions. The partial data of the total oxy- 
gen-isotope shifts are taken from Ref. jlTj and have been cor- 
rected for 100% oxygen-isotope exchange. The other data are 
from Refs. HH. 



FIG. 5. The site-selective oxygen-isotope effect in an un- 
derdoped Y0.7Pr0.3Ba2Cu3O6.97- The results show that the 
planar oxygen mainly (> 80%) contributes to the total oxy- 
gen-isotope shift in this underdoped cuprate. After Q. 

show the results of the site-selective oxygen-isotope effect 
for optimally doped YBa2Cu3 06.94. It is striking that 
the results shown in Fig. |] are in excellent agreement 
with those in Fig |[ This indicates a good reproducibil- 
ity of these experiments. The results for underdoped 
Y0.7Pr0.3Ba2Cu3O6.97 are plotted in Fig. M. Remark- 
ably, the planar oxygen ions mainly (> 80%) contribute 
to the total OIS in the underdoped samples as well. In 
Fig. ||, we show the total isotope shifts as a function of T c 
for Yi_ 2; Pr 2 ;Ba2Cu307_y together with the isotope shifts 
from the planar oxygen ions as well as from the apical and 
chain oxygen ions. From this figure, one can clearly see 
that the planar oxygen ions make a predominant contri- 
bution to the total OIS in all the doping levels. 



III. NEGLIGIBLE OXYGEN-ISOTOPE EFFECT 
ON THE CARRIER DENSITY 

From the above results, an important question arises: 
Are the observed large isotope shifts in underdoped 
cuprates caused by a possible difference in the carrier 
densities of two isotope samples or by a strong electron- 



There are three indirect experiments which have 
demonstrated that the difference in the hole densities of 
the 16 and 18 samples is smaller than 0.0002 per Cu 
site |p5[-p7f . Now we are able to determine the oxygen 
content very accurately using a very precise volumetric 
analysis |37) . Fig. (?] shows the oxygen contents of the 16 
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FIG. 7. The oxygen contents of the le O and ls O samples 
of La2-xSr :!; CuOi, with different doping x (unpublished data). 
The oxygen content was determined from a very precise vol- 
umetric analysis |37| . The weight of each sample used for the 
analysis is about 500 mg. The oxygen content can be deter- 
mined more precisely for a higher x and a heavier sample. The 
oxygen contents of two isotope samples are the same within 
±0.0002 per Cu site. 
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and 18 samples of I^-aSr^CuC^ with different doping 
x. It is remarkable that the oxygen contents of two iso- 
tope samples are the same within ±0.0002 per Cu site. 
These experiments consistently show that the difference 
in the hole densities of the 16 O and 18 O samples is negligi- 
ble, so that the observed large isotope effects are intrinsic 
and caused by a strong electron-phonon interaction. 



IV. LARGE OXYGEN-ISOTOPE EFFECT ON 
THE EFFECTIVE SUPERCARRIER MASS 

Several groups (2l]j2^] noticed that there is an isotope 
effect on the diamagnetic Meissner signal in decoupled 
fine-grained samples. Zhao et al. |25|] carefully studied 
this effect and interpreted it as due to the isotope-mass 
dependence of the average supercarrier mass m** (= 
V ( m ab) 2m ** )■ Since the magnetic penetration depth 
A(0) is proportional to y/m**/n s , then 

Am**/m** = 2AA(0)/A(0) + An s /n s , (3) 

where A means any small change of a quantity upon 
isotope substitution. Thus the isotope dependence of m** 
can be determined if one can independently measure the 
isotope dependence of A(0) and of n s . As shown above, 
there is a negligible oxygen-isotope effect on the normal 
carrier density n, and n s = n for clean superconductors, 
so one should expect no significant isotope effect on n s . 
Indeed, there is a negligible oxygen-isotope effect on the 
supercarrier density n s in YBa2Cu306.94 [ p8| . 

The isotope dependence of A(0) can be determined 
from that of the Meissner fraction /(0) which, for de- 
coupled and fine-grained samples, depends on the pen- 
etration depth A(0) and on the average grain radius R, 
as seen from the Shoenberg formula for spherical grains 



f(T) = -[l 



3(^)coth( 



R 
Mf) 



)+3( 



(4) 



where A(T) = \/[(A a t(T)] 2 A c (T) for layered compounds 
p0| . From equation EL it is obvious that a change in 
A(0) will lead to a change in /(0), which implies that the 
isotope dependence of A(0) can be determined from the 
isotope dependence of /(0). 

In Fig. ||, we show the Meissner effects for the 16 O and 
18 samples of La2- x Sr x Cu04 with x = 0.06 and 0.105. 
The samples are loosely packed with rather small grain 
sizes (R ~ 2-4 fim). One can clearly see that there are 
large oxygen isotope effects on both T c and the Meissner 
fraction. The most remarkable result is that the Meissner 
fraction of the 18 O sample is lower than that for the 16 O 
sample by about 23% in the case of x = 0.06. The iso- 
tope effects are reversible upon the isotope back-exchange 
[E6LE7|, and reproducible in several sets of samples. For x 



= 0.06, we evaluate Am**/m** = 24(2)%. This indicates 
a large negative oxygen-isotope effect on the effective su- 
percarrier mass in the deeply underdoped cuprates. 



5.0 
0.0 
5- -5.0 

p -10.0 

-15.0 

-20.0 

-25.0 




: 0.105 



• 16 

o 1 8q 



AT =-1.5(1) K 
Af/f = -6(1)% 



0000000° .• 
«••••••• 




10 15 20 25 30 



I I 

X = 


1 / '• ••••• ' 

0.06 Ap%.»»* 
/ i 


_ • 16 








o°°7 


J AT c = - 1.00(5) K 








Af/f = - 23(2)% 
i i i 



2 4 6 8 10 

Temperature (K) 

FIG. 8. The Meissner fractions for the 16 and ls O samples 
of La 2 - I Sr :l; Cu04 with x = 0.06 and 0.105. After [EiR. 



The above results were obtained from powder sam- 
ples and correspondingly reflect the average properties 
of those highly anisotropic superconductors. To gain 
more insight, it is essential to determine the oxygen- 
isotope effect on the inplane effective supercarrier mass 
m* I . Although we have tentatively extracted the isotope 
dependence of m*£ from experiments on powder samples 
p6[ , more reliable results should be obtained from exper- 
iments on single crystals. 

Unfortunately, a complete oxygen-isotope exchange by 
gas diffusion is impossible in single crystals with a large 
volume. In order to reach a complete oxygen-isotope ex- 
change, microcrystals with a volume of V ~ 150 x 150 x 
50 /tzm 3 (mass rj 10 /xg) should be used. In this case, 
commercial SQUID magnetometers do not have enough 
sensitivity to measure the magnetization for such tiny 
crystals especially near T c . Fortunately our highly sensi- 
tive torque magnetometer jllj is able to detect the small 
diamagnetic signal for the tiny crystals. 

The superconducting transition was studied by cooling 
the sample in a magnetic field B a — 0.1 T applied at an 
angle of 45° with respect to the c-axis 0]. The torque 
signal was continuously recorded upon cooling the sample 
at a cooling rate of 0.01 K/s. In order to determine the 
background signal of the cantilever, the measurement was 
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repeated in zero field and the data were subtracted from 
those of the field cooled measurement. The magnetic 
torque vs temperature data obtained for crystals with x 
= 0.086 and 0.080 are shown in Fig. ||. It is clear that T c is 
lower for the 18 exchanged samples. Furthermore, the 
magnetic signals of the back-exchanged samples (cross 
symbols) coincide with those of the 16 O annealed samples 
(open circles). We define T c as the temperature where 
the linearly extrapolated transition slope intersects the 
base line. The relative changes in T c are found to be 
AT C /T C = -5.5(4)% for x = 0.080 and -5.1(3)% for x 
= 0.086. The exponent ao is found to be 0.47(2) for 
x = 0.080 and 0.40(2) for x = 0.086, which is in good 
agreement with the results obtained for powder samples 
with similar doping levels |20|j2^]. 
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FIG. 9. The magnetic torque as a function of temperature 
for the le O and ls O microcrystals of La2- a; Sr :E Cu04 with x 
= 0.086 and 0.080. The measurements were carried out in a 
magnetic field B„ = 0.1 T applied at an angle of 45° with 
respect to the c-axis. After p2[. 

More interestingly, the inplane penetration depth 
A a b(T) can be extracted from field-dependent measure- 
ments Q. Fig. 10 displays A~ fc 2 (T) for the isotope- 
exchanged crystals with x — 0.086 and 0.080. The 
temperature dependence is described by the power law 
Kb( T ) = ^(O)! 1 - (T/T c ) n ] with an exponent n«5. 
From Fig. |icj it is evident that both T c and A~ h 2 (0) 
shift down upon replacing 16 O by 18 O. The shifts are 
found to be A\J{0)/\J(0) = -9(3)% and -7(1)% for 
x = 0.080 and 0.086, respectively. Using An s = 0, we 
find Am* a * b /m* a * b = 9(3)% for x = 0.080, and 7(1)% for x 
= 0.086. 
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FIG. 10. The inplane penetration depth A~ b 2 (T) for the 
16 and ls O microcrystals of L^-iSr^CuO^ with x = 0.086 
and 0.080. The A~ 6 2 (T) value was extracted from the field 
dependence of the torque. After 

Moreover, a substantial value of IS.m** b /m** b = 5- 
6% was also obtained for several optimally doped 
cuprates (e.g., YBa2Cu30g.94, Lai.ssSro.isCuC^, and 
Bii.6Pbo.4Sr2Ca2Cu30io+i;) from magnetization mea- 
surements p3| . Very recent muon-spin rotation exper- 
iments on the oxygen-isotope exchanged YBa2Cri306.96 
have confirmed that Am*l/m^l ~ 5%. Therefore, there 
is a substantial isotope effect on the in-plane supercarrier 
mass in optimally doped cuprates, but a very small ef- 
fect on T c . Such unusual isotope effects place strong con- 
straints on the pairing mechanism of high-temperature 
superconductivity. 



V. HUGE OXYGEN-ISOTOPE EFFECT ON THE 
CHARGE STRIPE FORMATION 
TEMPERATURE 

One of the most remarkable findings in the high- 
temperature copper oxide superconductors is the forma- 
tion of alternating spin and charge stripes below a char- 
acteristic temperature Q,^5|. Various x-ray absorption 
spectroscopic measurements 46 suggest that the local 
structures in the alternating stripes are different, form- 
ing an incommensurate superlattice. Such a stripe phase 
is believed to be important to the understanding of the 
pairing mechanism of high-temperature superconductiv- 
ity p7| . However, the microscopic origin of the stripe 
phase is still highly debated. It could be caused by purely 
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electronic interactions []47| and/or by a strong electron- 
phonon interaction 

Although there is increasing experimental evidence 
for a strong electron-phonon interaction in the cuprate 
superconductors, it is not clear whether this interac- 
tion is important to the formation of the stripe phase. 
For the colossal magnetoresistive manganites, a strong 
electron-phonon interaction plays an essential role in 
the formation of the Jahn- Teller stripes (or charge or- 
dering), as inferred from a very large oxygen-isotope 
shift of the charge-ordering temperature observed in 
both Ndo.5Sr .5Mn03 and Lao.sCao.sMnOa systems p9[ . 
Therefore it seems natural to seek for an isotope effect 
on the stripe formation temperature in cuprates. 
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techniques [Q. In Fig. |Tl|, we show the temperature 
dependence of the parameter R for the oxygen isotope 
exchanged Lai.g4Sro.o6Cu04. From the figure, one can 
clearly see that, upon replacing 16 with 18 0, the charge- 
stripe formation temperature T* in this cuprate increases 
from about 110 K to 180 K. Such a large negative iso- 
tope effect on charge ordering can be explained in terms 
of electron-phonon coupling beyond the Migdal approxi- 
mation [ pT| . 

VI. OXYGEN-ISOTOPE EFFECTS ON THE 
PSEUDOGAP FORMATION TEMPERATURE 




Temperature (K) 



FIG. 11. The temperature dependence of the parameter R 
for the oxygen-isotope exchanged Lai.g4Sro.o6Cu04. See text 
for the definition of the parameter R. After [ po[ . 

X-ray absorption near edge spectroscopy (XANES) is 
a powerful technique to probe the local structure con- 
formations for a system. From the XANES data, two 
characteristic peaks denoted by A and B in the XANES 
spectra of the cuprates are identified, which characterize 
the local structures within and out of the CuC>2 planes. 
A parameter R is defined as R — ([3i — ai)/(/?i + ai), 
where f3\ and a\ are the intensities of peak A and peak 
B in the XANES spectra, respectively. When charge- 
stripe ordering occurs, a change in the local structure 
takes place, leading to a sudden increase in R below the 
charge-stripe formation temperature T* . The identifica- 
tion of T* by XANES has been tested in a compound 
Lai.875Bao.i25Cu04 where T* was determined by other 



FIG. 12. Temperature dependence of the intrinsic 
linewidth T(T) of the -> I^ 1 ' transition of Ho 3+ for the 
16 and ls O samples of HoBa 2 Cu 4 08. After Ref. ||. 

Recently, a neutron spectroscopic investigation of the 
isotope effect on the relaxation rate of crystal field ex- 
citations of Ho 3+ in HoBa2Cu40g has been carried out 
p2| . In Fig. |l2], we plot the temperature dependence of 
the intrinsic linewidth T(T) (corresponding to the — > 
r^ 1 -* transition of Ho 3+ ) for the 16 and 18 samples. In 
the high temperature range, the linewidth appears to ex- 
hibit a linear temperature dependence. Cooling down to 
a characterisitc temperature T p , the linewidth suddenly 
gets narrower. Of particular interest is the characterisitc 
temperature T p which strongly depends on the oxygen- 
isotope mass. Upon replacing 16 by 18 0, T p changes 
from 170 K to 220 K, i.e., there is a large oxygen-isotope 
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effect on T p (« 50 K). Similarly, a large copper- isotope 
effect on T p has been found in the same compound |53] . 
The huge oxygen- isotope effects on T* 50 and on T p 
p2| indicate that a strong electron-phonon interaction 
plays an essential role in the charge dynamics and su- 
perconductivity |Q. Although the authors of Ref. |||] 
attributed the characterisitc temperature T p to a pseu- 
dogap formation temperature, one cannot rule out the 
possibility that T p may be related to a dynamical charge 
ordering. We would like to mention that a small positive 
oxygen-isotope effect on the spin pseudogap has been ob- 
served in YBa 2 Cu 4 8 pi. 



VII. OXYGEN-ISOTOPE EFFECT ON THE EPR 
LINEWIDTH 



O 



X 

< 




50 100 150 200 250 300 
T(K) 



3.0 



2.5 - 



O 

d. 2.0 



X 

< 



1.5 



1.0 











■ l8 o 


b) " 




□ l6 o 










a \" ■ ■ a 

h 

3 


□ 5.-°-' 


□ 



50 



100 



150 200 
T(K) 



250 300 



FIG. 



13. 



(a) 



Temperature depen- 
16 O isotope samples 



dence of the EPR linewidth AH PP for 
of La2-zSrzCuo.98Mno.02O4 with x = 0.03, 0.06, and 0.10; 
(b) The oxygen-isotope effect on the EPR linewidth AH PP in 
La1.g7Sro.03Cuo.98Mno.02O4. After Ref. |5^ , |59[ . 

The Electron Paramagnetic Resonance (EPR) is a 
powerful technique to probe and understand various phe- 
nomena in solid state physics. Unfortunately, the in- 
trinsic Cu 2+ signals in the cuprate superconductors have 
not been observed except for those of chain fragments 
in YBa 2 Cu 3 6+ 5 (0.7 < S < 0.9) and three-spin 



polarons |36|. The EPR silence has been attributed by 
many researchers to the very fast relaxation time of the 
Cu 2+ spins. In order to test this possibility, one needs to 
dope the material with a paramagnetic S-state ion (e.g., 
Mn 2+ ). Its relaxation occurs via the fast relaxing car- 
riers. This is called bottleneck effect, which has indeed 
been observed in Mn 2+ doped La2- a; Sr a ;Cu04 |57[] . 

Here we focus on the oxygen-isotope effect on the 
EPR linewidth in the Mn 2+ doped La2- a; Sr : j : Cu04 
p8| , p9[ . In Fig. [l^, we show the temperature depen- 
dence of the EPR linewidth AH pp for 16 O isotope sam- 
ples of La2-a;Sr 2; Cuo.98Mno.o204 with x = 0.03, 0.06, 
and 0.10 (Fig. |l5|a), as well as the isotope effect on 
the EPR linewidth AH pp in La1.97Sro.03Cuo.9sMno.02O4 
(Fig. [L3|b). It is remarkable that the EPR linewidth 
at 100 K for the 16 O sample is about half the one for 
the 18 O sample. Such a large oxygen-isotope effect on 
the linewidth indicates that the dynamics of the oxygen 
atoms play a key role in the Cu 2+ relaxation. Interest- 
ingly, the temperature dependencies of the linewidth for 
both isotope samples are in good agreement with th co- 
retical calulations (see solid and dashed lines) 
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VIII. LARGE OXYGEN-ISOTOPE EFFECT ON 
SPIN-GLASS FREEZING TEMPERATURE 
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FIG. 14. Temperature dependence of the internal magnetic 
fields (probed by muon-spin- rotation) for the 16 O and 18 O 
samples of La1.97Sro.03Cuo.9sMno.02O4. After 

The parent compounds of the cuprate superconduc- 
tors exhibit long-range 3D antiferromagnetic (AF) or- 
der, which is rapidly destroyed as holes are doped into 
the Cu02 planes. A short range ordered AF state exists 
at doping of 0.02<x<0.06 in La2- a; Sr :I ;Cu04. Early ^iSR 
and neutron scattering experiments found that this mag- 
netic state resembles a spin glass |Q. More detailed 
studies using 139 La nuclear quadrupole resonance ]6l| ] 
showed that the magnetic state in this doping regime 
is not a conventional spin glass, but a cluster spin glass 
(CSG). The understanding of how the short-range or- 
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dered AF state affects superconductivity and how it is 
influenced by lattice vibrations will help to clarify the 
pairing mechanism of high-T c superconductivity. 

It is well known that conventional theories of mag- 
netism neglect atomic vibrations; the atoms are mostly 
considered as being infinitely heavy in theoretical de- 
scriptions of magnetic phenomena, which rules out the 
existence of an isotope effect on magnetism. However, 
there is an exceptional example. In the ferromagnetic 
manganites, a giant oxygen-isotope effect on the Curie 
temperature has been observed |6^]. The question is 
whether such an isotope effect also exists in cuprates. 

In Fig. [l4|, we show the temperature dependence 
of the internal magnetic fields B p (probed by muon- 
spin rotation) for the 16 O and 18 O samples of 
La1.g7Sro.03Cuo.98Mno.02O4. It is remarkable that the 
spin-glass freezing temperature T g almost doubles upon 
replacing 16 O by 18 O. Such a huge isotope effect on the 
spin-glass freezing temperature suggests that spin dy- 
namics in cuprates is ultimately correlated with lattice 
vibrations. 



IX. OXYGEN-ISOTOPE EFFECT ON THE AF 
ORDERING TEMPERATURE 

The antiferromagnetic order observed in the parent 
insulating compounds like La2Cu04 signals a strong 
electron-electron Coulomb correlation. On the other 
hand, the large isotope effects found in the undcrdoped 
cuprate superconductors indicate a strong electron- 
phonon interaction. Now a question arises: can the 
strong electron-phonon interaction modify the antifer- 
romagnetic exchange energy and thus the AF ordering 
temperature in the parent insulating compounds? Stud- 
ies of the isotope effect on the AF ordering temperature 
could clarify this issue p4| , |65| ]. 

Fig. [l^ shows the temperature dependence of the 
susceptibility for the 16 O and 18 O samples of un- 
doped La2Cu04 (upper panel), and of oxygen doped 
La2Cu04+j, (lower panel). One can see that the AF or- 
dering temperature Tjv for the 18 O sample is lower than 
the 16 sample by about 1.9 K in the case of the un- 
doped samples. For the oxygen-doped samples, there is 
a negligible isotope effect. 

It is known that the antiferromagnetic properties of 
La2Cu04+j, can be well understood within mean- field 
theory which leads to a Tjv formula Bq] : 



k B T N = j'[£(T N )/a] 



(5) 



where J' is the inter layer coupling energy, £(Tjv) is the 
in-plane AF correlation length at Tjv with £(Tjv) oc 
exp(J/T/v) for y = ( J is the in-plane exchange energy). 
When Tjv is reduced to about 250 K by oxygen doping, 
a mesoscopic phase separation has taken place so that 



£(Tjv) = L (Ref. |37|), where L is the size of the anti- 
ferromagnetically correlated clusters, and depends only 
on the extra oxygen content y. In this case, we have 
Tjv = J'(L/a) 2 . Since L is independent of the isotope 



mass, a negligible isotope shift of Tjv in the oxygen-doped 
La2Cu04 +y suggests that J' is independent of the isotope 
mass. Then we easily find for undoped compounds 



AT N /T N = (AJ/J) 



B 



1 + B' 



(6) 



where B = 2T/T/v ~ 10. From the measured isotope 
shift of Tjv for the undoped samples, we obtain AJ/J 
~ -0.6%. 
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FIG. 15. The temperature dependence of the susceptibility 
for the 16 and ls O samples of undoped La2Cu04 (upper 
panel), and of the oxygen-doped La2Cu04+ H (lower panel). 
After Q. 

Theoretically, Kugel and Khomski |}8| considered the 
Jahn- Teller effect in a single-band Hubard model. They 
showed that, for U >> E p , Ticu, 



2f 2 



2E p hu> 
IP 



On the other hand, when U << Flu, 



J 



2t 2 exp(-2E p /Kuj) 
U -2E V ■ 



(7) 



(8) 



Here U is the onsite Coulomb repulsion, t is the bare 
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hopping integral, E p is the Jahn- Teller stabilization en- 
ergy, and ui is the vibration frequency of the Jahn- Teller 
mode. Assuming t = 0.5 eV, J = 0.13 eV, we obtain U 
= 3.8 eV. With hw = 0.1 eV, E p = 1.2 eV §§, we obtain 
A J I J ~ -0.1% from Eq. g which is about a factor of 6 
smaller than the measured value. The discrepancy might 
be due to the fact that the single-band Hubard model is 
oversimplified for the parent cuprates. 

X. CONCLUDING REMARKS 

In summary, the unconventional isotope effects ob- 
served in cuprates clearly demonstrate that the electron- 
phonon interaction plays an important role in the physics 
of cuprates. Our results also show that the phonon modes 
related to both oxygen and copper vibrations are impor- 
tant to the pairing. 

Interestingly, the concept of enhancing the electron- 
phonon coupling was the original motivation for the high- 
Tc discovery. Bednorz and Muller Q argued that there 
should be a strong electron-phonon interaction in per- 
ovskites with strong Jahn- Teller centers. Indeed, the 
stretching vibration mode, which is related to the Q2- 
type Jahn- Teller distortion, has proved to couple strongly 
to the doped holes |7fj|-f72|. This high-energy Q2-type 
mode may couple to the low-energy tilting mode {Q4/Q5- 
like Jahn- Teller mode) if the Cu-O-Cu bonding angle is 
less than 180°. Thus both stretching and tilting modes 
are important to high-T c superconductivity. 

Acknowlegement : The authors would like to thank K. 
A. Muller, A. Shengelaya, J. Hofer, A. Bianconi and D. 
E. Morris for their collaborations and discussions. The 
work was supported by the Swiss National Science Foun- 
dation. 



[1] J. G. Bednorz and K. A. Muller, Z. Phys. B 64, 189 
(1986). 

[2] B. Batlogg et al, Phys. Rev. Lett. 58, 2333 (1987). 
[3] L. C. Bourne et al, Phys. Rev. Lett., 58, 2337 (1987). 
[4] Donald E. Morris et al., Phys. Rev. B 37, 5936 (1988). 
[5] J.R. Schrieffer, X.G. Wen, S.C. Zang, Phys. Rev. B 39, 
11663 (1989). 

[6] A.J. Millis, H. Monien, D. Pines, Phys. Rev. B 42, 167 
(1990). 

[7] P.W. Anderson, The Theory of Superconductivity in the 
High-T c Cuprate Superconductors, Princeton University 
Press, 1998. 

[8] J.P. Carbotte, Rev. Mod. Phys. 62, 1027 (1990). 



[9] J. G. Daunt and K. Mendelssohn, Proc. Roy. Soc. A 185, 

225 (1946). 
[10] E. Maxwell, Phys. Rev. 78, 477 (1950). 
[11] C. A. Reynolds, B. Serin, W. H. Wright, and L. B. Nes- 

bitt, Phys . Rev. 78, 487 (1950). 
[12] H. Frohlich, Phys. Rev. 79, 845 (1950). 
[13] L. N. Cooper, Phys. Rev. 104, 1189 (1956). 
[14] M. R. Schafroth, Phys. Rev. 100, 463 (1955). 
[15] J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. 

Rev. 108, 1175 (1957). 
[16] A. S. Alexandrov and N. F. Mott, Polarons and Bipo- 

larons (World Scientific, Singapore, 1995). 
[17] J. P. Franck, in: Physical Properties of High Temperature 

Superconductors IV (World Scientific, Singapore 1994), 

ed. by D. M. Ginsberg, p. 189. 
[18] K. A. Muller, Z. Phys. B: Condens. Matter 80, 193 

(1990). 

[19] G. M. Zhao, K.Conder, M.Angst, S. M. Kazakov, J. 

Karpinski, M. Maciejewski, C. Bougerol, J. S. Pshirkov, 

and E. V. Antipov, Phys. Rev. B 62, Rll 977 (2000). 
[20] M. K. Crawford, M. N. Kunchur, W. E. Farneth, E. M. 

McCaron III, and S. J. Poon, Phys. Rev. B 41, 282(1990). 
[21] H. J. Bornemann, D. E. Morris, H. B. Liu, and P. K. 

Narwankar, Physica C 191, 211 (1992). 
[22] J. P. Franck, S. Harker, and J. H. Brewer, Phys. Rev. 

Lett. 71, 283 (1993). 
[23] D. Zech, K. Conder, H. Keller, E. Kaldis, E. Liarokapis, 

N. Poulakis, and K.A. Muller, in Proceedings of the 

Int. Workshop on Anharmonic Properties of High-T c 

Cuprates, eds. D. Mihailovic, G. Ruani, E. Kaldis, and 

K. A. Muller (World Scientific, Singapore 1995) p. 18. 
[24] D. Zech, K. Conder, H. Keller, E. Kaldis, and K. 

A. Muller, Physica B 219 & 220, 136 (1996). 
[25] G. M. Zhao, K. K. Singh, A. P. B. Sinha, and D. E. 

Morris,, Phys. Rev. B 52, 6840 (1995). 
[26] G. M. Zhao, M. B. Hunt, H. Keller, and K. A. Muller, 

Nature (London) 385, 236 (1997). 
[27] G. M. Zhao, K. Conder, H. Keller, and K. A. Muller, J. 

Phys.: Condens. Matter, 10, 9055 (1998). 
[28] W. E. Pickett, R. E. Cohen, and H. Krakauer, Phys. Rev. 

Lett. 67, 228 (1991). 
[29] G. M. Zhao, V. Kirtikar, K. K. Singh, A. P. B. Sinha, D. 

E. Morris, and A. V. Inyushkin, Phys. Rev. B 54, 14956 

(1996). 

[30] D. E. Morris, A. P. B. Sinha, V. Kirtikar, K. K. Singh, 

and A. V. Inyushkin, Physica C 298, 203 (1998). 
[31] G. V. M. Williams, D. J. Pringle, and J. L. Tallon, Phys. 

Rev. B 61, R9257 (2000). 
[32] M. Cardona, R. Liu, C. Thomsen, W. Kress, E. 

Schonherr, M. Bauer, L. Genzel, and W. Konig, Solid. 

State Commun. 67, 789 (1988). 
[33] W. K. Ham, S. W. Keller, J. N. Michaels, A. M. Stacey, 

D. Krillov, D. T. Hodul, and R. H. Fleming, J. Mater. 

Res. 4, 504 (1989). 
[34] J. H. Nickel, D. E. Morris, and J. W. Ager III, Phys. Rev. 

Lett. 70, 81 (1993). 
[35] D. Zech, H. Keller, K. Conder, E. Kaldis, E. Liarokapis, 

N. Poulakis, and K. A Muller, Nature 371, 681 (1994). 
[36] G. M. Zhao, J. W. Ager III, and D. E. Morris, Phys. Rev. 

B 54, 14982 (1996). 



10 



[37] K. Conder, S. Rusiecki, and E. Kaldis, Mat. Res. Bull., 
24, 581 (1989). 

[38] G. M. Zhao and D. E. Morris, Phys. Rev. B 51, 16487 
(1995). 

[39] D. Shoenberg, Proc. R. Soc. London, Ser. A175, 49-70 
(1940). 

[40] V. G. Kogan, M. M. Fang, and S. Mitra, Phys. Rev. 

B38, 11958 (1988). 
[41] M. Willemin et al, J. Appl. Phys. 83, 1163 (1998). 
[42] J. Hofer, K. Conder, T. Sasagawa, G. M. Zhao, M. 

Willemin, H. Keller, and K. Kishio, Phys. Rev. Lett. 84, 

4192 (2000). 

[43] G. M. Zhao, V. Kirtikar, and D. E. Morris, Phys. Rev. B 
(in press). 

[44] J. M. Tranquada, B. J. Sternlleb, J. D. Axe, Y. Naka- 
mura, and S. Uchida, Nature (London) 375, 561 (1995). 
[45] H. A. Mook, P. C. Dai, S. M. Hayden, G. Aeppli, T. 

G. Perring, and F. Dogan, Nature (London) 395, 580 
(1998). 

[46] A. Bianconi, N. L. Saini, A. Lanzara, M. Missori, T. 
Rossetti, H. Oyanagi, H. Yamaguchi, K. Oka, and T. Ro, 
Phys. Rev. Lett. 76, 3412 (1996). 

[47] V. J. Emery, S. A. Kivelson, and O. Zachar, Phys. Rev. 
B 56, 6120 (1997). 

[48] A. Bussmann-Holder, A. Simon, H. Buttner, and A. R. 
Bishop, Phil. Mag. B 80, 1955 (2000); A. Bussmann- 
Holder, K. A. Midler, R. Micnas, H. Buttner, A. Simon, 
A. R. Bishop, and T. Egami, J. Phys.: Condens. Matter, 
13, L169 (2001). 

[49] G. M. Zhao, K. Ghosh, and R. L. Greene, J. Phys.: Con- 
dens. Matter, 10, L737 (1998); G. M. Zhao, K. Ghosh, 

H. Keller, and R. L. Greene, Phys. Rev. B 59, 81 (1999). 
[50] A. Lanzara, G. M. Zhao, N. L. Saini, A. Bianconi, K. 

Conder, H. Keller, and K. A. Miiller, J. Phys.: Condens. 
Matter 11, L541 (1999). 
[51] S. Blawid and A. J. Millis, Phys. Rev. B 63, 115114 
(2001). 

[52] D. R. Temprano, J. Mesot, S. Janssen, K. Conder, A. 

Furrer, H. Mutka, and K. A. Miiller, Phys. Rev. Lett. 

84, 1990 (2000). 
[53] D. R. Temprano, J. Mesot, S. Janssen, K. Conder, A. 

Furrer, A. Sokolov, V. Trounov, S. M. Kazakov, J. 

Karpinski, and K. A. Miiller, Eur. Phys. J. B 19, 5 

(2001). 

[54] F. Raffa, T. Ohno, M. Mali, J. Roos, D. Brinkmann, K. 
Conder, and M. Eremin, Phys. Rev. Lett, bf 81, 5912 
(1998). 

[55] J. Sichelschmitt et al, Phys. Rev. B 51, (1995) 

[56] B.I. Kochelaev et al, Phys. Rev, Lett. 79, 4274 (1997). 

[57] B.I. Kochelaev et al, Phys. Rev. B 49, 13106 (1994). 

[58] K. A. Miiller, Physica C 341-348, 11 (2000). 

[59] A. Shengelaya et al, Phys. Rev. B 63, 144513 (2001). 

[60] B. J. Sternlieb, G. M. Luke, Y. J. Uemura, T. M. Rise- 
man, J. H. Brewer, P. M. Gehring, K. Yamada, Y. Hi- 
daka, T. Murakami, T. R. Thurston, and R. J. Birgeneau, 
Phys. Rev. B 41, 8866 (1990). 

[61] J. H. Cho, F. Borsa, D. C. Johnston, and D. R. Torgeson, 
Phys. Rev. B 46, 3179 (1992). 

[62] G. M. Zhao, K. Conder, H. Keller, and K. A. Miiller, 
Nature (London) 381, 676 (1996). 



[63 

[64 
[65 



[66 
[67 

[68 

[69 
[70 



[71 
[72 



A. Shengelaya, G. M. Zhao, C. M. Aegerter, K. Conder, 
I. M. Savic, and H. Keller, Phys. Rev. Lett. 83, 5142 
(1999). 

G. M. Zhao, K. K. Singh, and D. E. Morris, Phys. Rev. 
B 50, 4112 (1994). 

S. Kramer, K. H. Gyorfi, A. Bussmann-Holder, K. Con- 
der, and M. Mehring, Phys. Stat. Solid. B 215, 601 
(1999). 

T. Thio et al., Phys. Rev. B 38, 905 (1988). 

J. H. Cho, F. C. Chou, and D. C. Johnston, Phys. Rev. 

Lett. 70, 222 (1993). 

K. I. Kugel and D. I. Khomski, Sov. Phys. JETP 52, 501 
(1980). 

H. Kamimura, Int. J. Mod. Phys. B 1, 873 (1987). 

L. Pintschovius and W. Reichardt, in Physical Properties 
of High Temperature Superconductors IV, edited by D. 
Ginsberg (World Scientific, Singapore, 1994), p. 295. 
R. J. McQueeney, Y. Petrov, T. Egami, M. Yethiraj, G. 
Shirane, and Y. Endoh, Phys. Rev. Lett. 82, 628 (1999). 
Y. Petrov, T. Egami, R. J. McQueeney, M Yethiraj, H. 
A. Mook, and F. Dogan, hond-mat/0003414. 



11 



